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ABSTRACT 

We present Spitzer Space Telescope and Herschel Space Observatory infrared observations 
of the recurrent nova T Pyx during its 201 1 eruption, complemented by ground-base optical- 
infrared photometry. We find that the eruption has heated dust in the pre-existing nebulosity 
associated with T Pyx. This is most likely interstellar dust swept up by T Pyx - either dur- 
ing previous eruptions or by a wind - rather than the accumulation of dust produced during 
eruptions. 

Key words: circumstellar matter - stars: individual, T Pyx - novae, cataclysmic variables - 
infrared: stars - ISM: general 



1 INTRODUCTION 

Nova eruptions occur as a result of a thermonuclear runaway (TNR) 
on the surface of a white dwarf (WD) in a semi-detached binary 
system. In classical novae (CNe) the secondary star is normally a 
red dwarf. Mass is transferred from the secondary through the in- 
ner Lagrangian point onto the surface of the WD via an accretion 
disc. The degenerate layer o f accreted material is compress ed and 
heated, and a TNR occurs dStarrfield. lUadis & HixllBosI) . Con- 
sequently ~ 10^'' Mq of material, enriched in CNO (and other 
metals), is ejected at a few hundred to a few thousand kms~^ 
jGehrz et al.lll998l : iBode & Evans! |2008|) in a CN eruption. Once 

* Herschel is an ESA space observatory with science instruments provided 
by European-led Principal Investigator consortia and with important partic- 
ipation from NASA, 
t E-mail: ae@astro.keele.ac.uk 



the eruption has subsided the mass transfer resumes, and in time 
(~ 10* — 10^ years) another nova eruption occurs. A CN system 
therefore undergoes many eruptions during its lifetime: nova erup- 
tions recur. 

Recurrent novae (RN) undergo the same evolution; however 
eruptions r ecur on less than a human timescale, typically 10-20 
years (e.g. lEvans et al.ll2008|, for a recen t review). They seem to 
divide into 3 sub-classes ( lAnupamal l2008'). These are the "U Sco" 
class (with short orbital periods a nd spectral evolu tion similar to 
that of the 'He/N' class of CNe; see IWilliamsll99"3) : the "RS Oph" 
class (with long orbital periods), and the "T Pyx" class (also 
with short orbital periods and spectral evolution that evolves from 
'He/N' to 'Fen'). The T Pyx class resembles the CNe in terms of 
the nature of the binary system and the spectral evolution during 
eruption. 

We present here infrared (IR) space- and ground-based IR 
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Figure 1. The visual (V-band) light curve of T Pyx during the 201 1 erup- 
tion, based on observations in the AAVSO (black) and SMARTS (orange) 
databases. The apparent offset for t 150 days arises because AAVSO data 
have been corrected for a differential airmass effect, caused by most of the 
flux being in str ong discrete lines rather tha n being smoothly distributed in 
the continuum iOksanen & SchaefeilEoT l|). The times of the IR observa- 
tions discussed in this paper are indicated by the arrows. 



spectrophotometry, as well as optical broadband photometry, of the 
RN T Pyx following its 20 11 eruption. 



2 TPyx 

2.1 The system 

TPyx has undergone recorded nova eruptions in 1890, 1902, 
1920, 1944 and 1966; in fact, it was the first nova identified as 
a recurrent. The T Pyx binary has an orbital period of 1.83 h 
dUthas. Knigge & St eeghs 20 id ). The mass ratio is determined to 



a distance > 4.5 kpc and reddening E{B - V) = 0.50 ± 0.10. 
We shall use D = 4 kpc and E{B — V) = 0.4 here, although the 
precise values do not affect our main conclusions (e.g., reddening 
has negligible effect at the IR wavelengths of interest here). 

2.2 The 2011 eruption 

T Pyx was discovered to be in ou tburst on 2011 April 14.29 UT 
(MJD55665.29) by M. Linnolt (see lWaageij|201 ll ). The 201 1 erup- 
tion, al most 50 years after its last, was well overdue; this de- 
lay led iSchaefer et al.l ( I2OIOI) to speculate that T Pyx was enter- 
ing "hibernation", and that there would be no further eruptions for 
~ 10" years . Possible reasons for this inter-outburst behaviour are 
discussed bv lSchaefer et al. 1( 1201 1), who also provide a comprehen- 
sive discussion of the early visual l ight curve. 

Spectroscopic observations bv lShore et al.l (1201 lb revealed ex- 
pansion velocities ^ 2 500 km s~^ during the early stages of the 
eruption. They deduce an ejected mass of 10~^ /M©, where / is 
the filling factor for the ej ected material. VLTI a nd Mount Wilson 
CHARA IR observations jChesneau et alj[201 ih provide evidence 
for a bi-polar ejection that is e ssentially face-on, consistent with the 
low inclination of the system fa thas et alj2010l) . 

Swift (Gehrels 2004) observations started som e 7.5 hours af- 
ter dis covery, and revealed a soft X-ray source (iKuuIkers et al.l 
I2OI lal) , but between days 12-109 (which c overs the IR data re - 
ported here) the X-ray source was very weak dOsbome et al.bOllh . 
Subsequently (from day 142), the X-ray emission was bright and 
variable (Kuulkers et al. 2011b). 

The visual (V^-band) light curve during the 2011 eruption is 
shown in Fig. [T] The UVOT instrument on Swift showed that the 
ultra violet flux at 2246A peaked some 25 days later than the visual 
flux jOsbomeetalltollh . 



be 0.20 ± 0.03 by Uthas et aljfeoioh , who estimate the WD and 3 OBSERVATIONS 



cool companion mass to be 0.7 M© and 0. 14M0 respectively. The 
former is rather low for either a RN or a CN and if the WD mass is a 
more plausible ^ 1 M0 , the secondary mass has to be ^ 0.2 M© ; in 
the discussion be l ow we assume 1.2 Mq for the WD mass in T Pyx 
( lAnupamj|200i ; ISchaefer et id]|2010l) . lUthas et al.l ( I2OIOI) deter- 
mine that the inclination of the binary is ~ 10°, i.e. nearly face- 
on. Possible e volutionary scenarios fo r T P yx have been discu ssed 
by Schaefer. Pagnotta & Shitfal ( 1201 Ol) and lUthas et al.l ( I2OIOI) . In 
T Pyx, as in a CN, the nova explosion occurs on the surface of the 
WD following a TNR in a layer of degenerate material accreted 
from the secondary. 

T Pyx is unusual amongst RNe i n that it lies at the centre of a 
nova-like shell with diameter ~ 15 dDurbeckfe Seitteilll979l see 
alsolSchaefer et al1l201(]t) . A detailed study of knots in the nebula 
( ISchaeferetal.112010! )" has revealed complex interactions between 
material ejected in earlier eruptions. Schaefer et al. find that there 
has been no significant deceleration of the knots and, based on the 
angular expansion rate of the knots, deduce that they were ejected 
by a CN-type eruption close to the year 1866. They find knots 
ejected in 1866 that have "turned on" since 1995, and suggest that 
these knots are powered by shocks caused by collisions with fast 
ejec ta from more r e cent R N eruptions. 

ISelvelli et aL bOOSi) est i mated the distance of T Pyx as 
3500 ± 350 pc. lSchaefer et all ( I2OI0I) summarise several estimates 
of the interstellar reddening, and conclude that E[B — V) = 
0.25 ± 0.02. On the basis of observati ons carried out in the after- 
math of the 201 1 eruption (see below), IShoreetal.l(l201lh propose 



Broadband photometric IR observations w ere carried out at 
the Infrared Telescope on M ount Abu dDeshpand^ 1 19951 : 
lAnandarao & Chakrabortvl I2OI0I) . and as Directors' Discre- 
tionary Time on t he Spitzer Space Telescope dWemer et aL 20041; 
Gehrz et l2007h and the Herschel Space Observatory dPilbratd 
20031 ; IPilbratt et al.l I2OI0I) ; these data were supplemented, as 



closely in time as possible, by optical and IR photometry from the 
AAVSCQ and SMART^J archives. 



3.1 Mt. Abu 

Near IR photometry of T Pyx in the JHK bands was carried out 
with the Mt. Abu 1.2-m telescope using the 256 x 256 HgCdTe 
NICMOS3 array of the Near- Infrared Imager/Spectrometer (see 
iBaneriee & Asho^l201lh . Observations of both T Pyx and a com- 
parison star (SAO 177754: J = 7.082±0.029, H = 7.070±0.061, 
K = 6.963 ±0.018) were made at 5 dithered positions in each of 
the J, H and K bands. The dithered images were combined to pro- 
duce a median sky frame, with dark counts included, which was 
subsequently subtracted from the object frames. 



^ http://www.aavso.org/ 

^ http://www.astro.sunysb.edu/fwalter/SMARTS/NovaAtlas/ 
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Table 1. Optical and infrared fluxes for T Pyx. Ground-based data are given both in magnitudes and, immediately beneath, in mJy; IRAC and PACS data are in mjy. All mJy values are in bold font. £ is the time (in 
days) since eruption. Numbers in brackets are uncertainties in last figure(s) of quoted fluxes. Horizontal lines separate quasi-contemporaneous data obtained on different dates. 
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PACS (mJy) 



YYYY-MM-DD.DD 


MJD 


t(d) 


Facility* 
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V 
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J 


H 


K 


3.6 ^m 


4.5 ^m 


70 fj.m 


100 ^m 


130 ^tm 


2011-05-09.52 


55690.52 


25.23 
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7.47(1) 


7.18(1) 


6.56(1) 


6.26(1) 


- 


- 


4.97(1) 
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- 


- 


- 
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- 
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- 
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- 


- 
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- 
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55716.51 


51.22 


1 
















2048(10) 


2000(10) 








2011-06-04.96 


55716.96 
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*A = AAVSO; B = Herschel FACS; I = Spitzer IRAC; M = Mt Abu; S = SMARTS. 
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Aperture photometry of the sky-subtracted frames was done using 
IRAF to yield the JHK magnitudes given in Table[T] the magnitudes were 
converted to Jy using standard zero magnitude fluxes in ICoxl teOOOi) . De- 
tailed results of the near-IR studies of T Pyx from Mt. Abu wiU be presented 
elsewhere. 



3.2 Spitzer Space Telescope 

T Pyx was observed as a target of opportunity (ToO) target using a Di- 
rector's Discretionary Time allocation (PID: 70260) with the Infrared 
Array Camera (IRA C; iFazio et al] |2004|) on the Spitzer Space Telescope 
iWemer et al.ll2004l) at wavelengths 3.6 ^m and 4.5 pm, on two occasions 
as detailed in Table[T] The observations were carried out in full array mode, 
and consisted of 9-point random, medium scale dither patterns. The frame 
time was 0.4 s and the on-source integration time was 3.6 s per pointing. 
The data were reduced using MOPEX CMOPEX.201 li) . 

3.3 Herschel Space Observatory 

T Pyx was observed with the Photodet ector A rray Camera and 

Spectrometer instrument ( PACS; Poglitsch et al ] I2OIOI) on the 
Herschel Space Observatory jPilbratt et alj I2OIOI) . Photometry with 
PACS was obtained in Line Scan mode; the area covered was typically 
315" X 237" and the on-source integration time was 110 s. Spectroscopy 
was also obtained using PACS in Line Spectroscopy mode, to cover the 
lines Hi (11-10) 52.5 /im, [NlII]57.32^m, [OlII]51.80^m, Hi (15-13) 
61.9 /im, [Ol] 63.18/tm, [N IV] 69.40 ^m, [Ov]73.5 ^m, [N II] 76.50 ^tm, 
[O III] 88.36 /im and [N IV] 158.5 /.tm. Integration times ranged from 
2160 s to 3 072 s on-s ource. All the PACS data were reduced using HIPE 
version 7.0 l lOtll2010l) . 

Fluxes in the imaging observations were measured using an object 
aperture of 48" and the background was measured typically in an annulus 
of 96" (inner) and 144" (outer) diameter. T Pyx was securely deteced by 
PACS in Photometry mode; the fluxes are given in Table [T] Apart from 
[N III] 57.32 pm (which we will discuss elsewhere), none of the above hnes 
were detected in spectroscopy mode in the Herschel aperture, to 3cr limits 
of typically ~ 100 mJy in all three bands (see Fig. |2). The [NlII] line is 
at the edge of the 70 /im band and will not contribute significantly to the 
observed flux. 

We note that T Pyx is unresolved in both IRAC and PACS images. All 
the data are summarised in Table[T] 



4 RESULTS AND DISCUSSION 

For wavelengths ^2.2 pm the spectral energy distiibution (SED) of 
T Pyx resembles a black body with temperature 9 000 K. There is 
a clear excess at long wavelengths (see Fig. |3] that appears in the 
Herschel Space Observatory PACS data. While there are several fine struc- 
ture lines that may contribute to the broadband PACS fluxes, our PACS 
spectroscopic observations show that emission lines can not have con- 
tributed to the observed PACS fluxes. 

We show below that the far-IR SED is consistent with emission at 
~ 45 K and we conclude that the far-IR emission we see is due to cool 
dust. Newly formed dust would have a te mperature o f several hundreds of 
degrees K so soon after the eruption (e.g.. lEvans et al . 2005), and would be 
evident in near-IR ( ^ 5 pm) spectra obtained throughout the early phase of 
the eruption (Woodward et al., in preparation). Furthermore, moving at ~ 
2 500 km s~^, the 201 1 ejecta would not reach a distance from the central 
star commensurate with the low dust temperature for ^ 100 years after 
outburst. This rales out the formation of the cool dust in the material ejected 
in the 201 1 eraption. 

The visual light curves of RNe tend to be replicated from outburst to 
outburst ISchaefer 2010); TPyx is no exception, and there is no evidence 
for extinction events in the fight curve that would be associated with dust 
formation during the eruptions of T Pyx. The lack of dust formation is very 
likely in the case of RNe like RS Oph (in which the ejected material runs 
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Figure 3. Spectral energy distribution of T Pyx; data are dereddened by 
E{B — V) = 0.4. The sofid and broken black curves are DUSTY fits, with 
Sil and AmC grains respectively, to the BVRI, JHK and PACS data for 
MJD 55691.5 (2011 May 10.53 UT; filled circles). The other curves are 
9 000 K blackbody fits to BVRI and IRAC data for MJD 55690.5 (201 1 
May 9.52; triangles) and MJD 55716.96 (201 1 June 4.96 UT; squares). See 
text for details. 

into and shocks the stellar wind; lEvans et"al]|20073) . However we take a 
more cautious approach in the case of T Pyx, particularly as the T Pyx class 
of RNe resemble CNe, many of which are known dust-producers. We there- 
fore do not rule out the possibility that dust might have formed in the ejecta 
of previous RN eruptions (even if it did not do so in the 201 1 eruption). 

We explore the origin of the dust evident in emission at A ^ 30 pm 
(Fig. |3) in terms of (i) dust ejected in previous T Pyx eraptions and (ii) in- 
terstellar dus t swept up by ejecta, or by wind s originating in the T Pyx 
system (e.g. iKnigge, King & Patterso j l200(i ). We consider two cases: 
(a) amorphous carbon (AmC) grains of radius a = 0.2 ^tm, such as con- 
dense in CN wind s (see e.g. lEvans et"ai]ll997l I2O0I lOehrz et al.lll998l ; 
lEvans & Rawlings 2008) and (b) 0.1 ^m silicate (Sil) grains. In case (a) 
the grains represent dust formed in previous eruptions and accumulated in 
the environment of T Pyx, while case (b) represents swept-up interstellar 

dust. 

We have modelled the SED using the DUSTY code jivezic & Elitzud 

I1995I) . using a 9 OOP K blac kbody at the Eddington luminosity for a 1 .2 M0 
WD )Schaefer et a"i]|2O10l) as the input s ource (see above); w e note that 
the X-ray emission at this time was weak jOsbome et al.ll201ll) and we are 
therefore justified in not including any other sources (especially hot sources) 
of dust heating in our modelling. We assume a geomet rically thin dust shell 
and take opt i cal co nstants for Sil grains from Draine & Led ( Il984h . and 
from lHannej h988h for AmC. 

We find that, in both AmC and Sil cases the dust temperature at the 
inner boundary of the dust shell is 45 it 5 K, signific antly higher than that of 
dust in the interstellar medium (~ 20 K: ICoxl200 d); we conclude therefore 
that the dust is indeed associated with T Pyx. We further find that the SED is 
well fitted by an optically thin dust shell (optical depth ty — 1.0[it0.2] X 
10-* at V for AmC, ry ~ 2.2[±0.3] X 10""' for Sil). The dust mass 
is ~ 8.3 X 10-5 M0 (AmC) or ~ 2.1 X 10-^ Mq (Sil); see Fig.[3]and 
Table|2]for other parameters). For a thin shell with inner radius ri the dust 
mass scales as r^r, and the value of ri is fixed by the assumed luminosity 
of the central source, the dust temperature at the inner boundary and the 
grain material. The main source of uncertainty is in r, and so the uncertainty 
in the dust mass (SM;just is given by <5M(just /A/dust — <5t/t. Table |2] 
summarizes the best-fit DUSTY model parameters. 

While the optical depth and temperature of the dust shell are reason- 
ably well constrained by the SED, the data are inadequate to provide any 
information about grain size or composition. However the deduced angular 
diameter of the dust shell is, for the Sil case, ~ 33", comparable with 
the dimensions of the optical nebulosity associated with T Pyx (and less 
than the aperture used to measure the Herschel fluxes; see above). The cor- 
responding value for AmC is substantially greater than this, ~ 78", and 
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Table 2. Parameters for DUSTY fits for a 1.2M0 WD; 9 is tlie angular diameter corresponding to the diameter 2ri of the dust shell. In both cases a dust shell 
with thickness 0.2ri has been assumed. See text for details and for explanation of scaled mass. 



Dust type 


T(K) 


a ( fim) 


ri (m) 


8 (arcsec) 


TV 


Dust mass ( Mq ) Scaled mass ( Mq ) 


AmC 


45 


0.2 


2.3 X IQiS 


78 


1.0 X 10-* 


8.3 X 10"^ 


Sil 


45 


0.1 


1.0 X lO^S 


33 


2.2 X 10"* 


2.1 X 10-^ 6.2 X 10-4 









e 


{ * 

1 








Observer 



The mass of interstellar dust contained in a sphere having the angular 
diameter of the d usty T Pyx ne bula (~ 33"), assuming 10~^ 0.1 fim sil- 
icate grains m-^ iAllenlll973[) . is expected to be ^ 3.1 X 10"^ Mq (the 
corresponding mass in the Herschel aperture is ~ 9.2 X 10^^ Mq). This 
is of the same order as the dust mass we detect. The simplest interpreta- 
tion therefore is that we are seeing the result of interstellar dust swept up 
by winds from T Pyx. As the dust shell will continue to be heated by the 
2011 emption for some time (> 2t\/c ~ 2 years) there is ample time 
to plan further far-IR observations to place tighter constraints on the shell 
parameters deduced here. 



Figure 4. Geometry of IR echo. Star at O "switches on" at time t = and 
illuminates dust shell (shaded); locus of constant light travel time at time t 
after eruption is parabola (dotted line) with focus at O and semi-latus rectum 
OS = ct/2. From the point of view of a distant observer only the hatched 
area is perceived to be illuminated; the observed diameter is labelled by d. 

significantly greater than the dimensions of the nebulosity; we restrict the 
remainder of the discussion to the Sil case. 

We should inject a cautionary note in that the inner radius of the dust 
shell as deduced from the DUSTY modelling is r\ ~ 1.0 X 10*'' m, or ~ 
390 light days. In view of the fac t that T Pyx is highly variable, light-travel 
times are likely to be important teode & Evansll 19791 see also Fig. |4) but 
these are not built into DUSTY. It can be shown that the angular diameter of 
the heated region as seen by an infinitely distant observer, att = 26.24 days 
(the time of the PACS observation) and for the above r\, is 

2ct I ( 2ri \" 

this is somewhat less than the "full" angular diameter of the dust shell (r^ 
33") and of the optical nebulosity (~ 15"). 

Also, if the (geometrically thin) dust shell is spheiically symmetric it 
can be shown that only a fraction ct/2ri of the dust shell is perceived by 
the observer to be illuminated (i.e., heated); at t = 26.24 days this fraction 
is 3.4% for Sil grains. DUSTY models of the cool dust emission yield a Sil 
dust mass of 2.1 X 10-^5 Mq (Table|2); however, only 3.4% of the total 
dust mass is actually perceived to be illuminated at this epoch. Thus a lower 
limit to the Sil total dust mass can be determined simply by scaling. We 
estimate that there is actually ~ 6.2 X 10^* Mq of silicate dust present. 

We can definitely rule out dust-fonnation in the 201 1 ejecta of T Pyx 
because there is no near-lR evidence. In addition, the "scaled mass" of dust 
far exceeds that produced in a nova eruption. Indeed, if confined to a cir- 
cumstellar shell of radius ~ 5.7x lO*'^ m (the distance that ejecta moving at 
2 500 km s" * would travel in 26.24 days), it would have resulted in substan- 
tial extinction in the optical, and the corresponding ejected mass (assuming 
a gas-to-dust ratio of ~ 100) would be very large. We conclude that the 
cool dust we see in the environment of T Pyx is either (a) condensed dust 
accumulated from earlier eruptions, or (b) interstellar dust that has been 
swept up by material ejected by T Pyx in previous eruptions o r by a wind 
from t he underlying binary. If the ejected mass determined bv lShore et alj 
i201 ll) . IQ— ^ /Mq, is typical of previous eruptions the dust-to-gas ratio is 
unfeasibly large. If the dust we see is dust condensed in previous eruptions 
it must be the accumulation of at least ~ 60 eruptions (as / ^ 1, and not 
all eruptions will have produced dust). 



5 CONCLUSIONS 

We have presented IR observations of the RN T Pyx following its 201 1 out- 
burst. We see a cool, weak IR excess, which we argue is due to the heating 
of dust in the environment of T Pyx that was present before the 201 1 erup- 
tion. We attribute this to interstellar dust that has been swept up by material 
ejected by the T Pyx system, either in the course of eruptions or a wind from 
the binary. 

It is particularly interesting that two RNe have now been found to 
have dusty environments that pre-date their eruptions, although the respec- 
tive circumstances are very different. While the dust around T Pyx is cool, 
lEvans et alj j2007b) found hot silicate dust in the immediate environment 
of RS Oph shortly after its 2006 outburst, and in this case the dust seemed 
to have survived the extreme environment (hard radiation field and shock 
blast wave) of the 2006 eruption (see also Rushton et al, to be submitted). 
On the other hand there is no e vidence for any dust in the environmen t 
of the RN U Sco in the near-IR tevans et alJl200ll ; lBaneriee et alJl20ioh . 
These three RNe represent very different facets of the RN phenomenon and 
indeed, each is representative of the specific classes of RN l lAnupamal2008l) 
noted in Section 1 , although membership of each class is very small. 

Further observations of RNe, particularly at wavelengths 10 ftm, 
would help to throw light on the environments of RNe, and hence poten- 
tially on their recent history and evolution. 
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